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Abstract 

The X-34 Main Propulsion System (MPS) configuration 
includes the liquid oxygen (LOX) and rocket propellant 
#1 (RP-I) feedlines (Figure 1) The flow analyses of 
these feedlines were performed and documented in 
previous studies [l \ These analyses predicted a relatively 
low inlet distortion and nearly even flow split at the 
engine interface. The new design for these MPS 
feedlines has been recommended recently. The new 
configuration includes a tighter radius in the RP-1 
feedline and a neck-down section between the gimbals 
(Figure 2). Conversely, the LOX feedline is very similar 
to the previous design. There were concerns that this 
new RP-1 configuration might generate a greater flow 
distortion at the engine interface than the original 
design. To resolve this issue, a Computation Fluid 
Dynamics (CFD) analysis was conducted to determine 
the flow field in the new RP-1 feedlines. 


Introduction 

The flow field at the engine interface can have a large 
effect on pump performance and structural integrity. 
Distorted flow at the pump inlet can lead to significant 
decreases in suction capability and head rise as well as 
unacceptably high dynamic loads. Because of this CFD 
analyses have been performed for the LOX and RP-1 
feedlines. 

This paper describes the procedure to obtain the flow 
field in the RP-1 feedline and reports the predicted 
environments at the engine interface as well as the 
simplifications and assumptions embedded in the 
analyses. The flow analyses of the new LOX feedline 
configuration will be included in the full paper. 

Modeling 

A model was created and analyses were performed on 
the new design of the MPS feedlines. The new analyses 
included the pressure-compensated bellow, which for 
simplification was disregarded in previous analyses. 
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Figure 3 shows the configurations of both designs. The 
new design, shown in wireframe, is overlaid on top of 
the old design. Note that the model does not include the 
gimbal joints, fuel filters, and flow meters. The 
computational domain for the new design included 3 
blocks with a total of 239,463 grid points. The grids 
were generated using the code GRIDGEN m . The 
solutions were obtained using the code FDNS l2 \ a 
three-dimensional Navier-Stokes flow solver, with a k-e 
turbulence model. The analyses performed were single 
phase, incompressible, and steady state. Fully developed 
velocity and turbulence profiles were used at the inlet 
and mass conservation was enforced at the exit. At the 
pipe surfaces, no-slip wall boundary conditions were 
applied. The flow fields were initialized to the reference 
conditions, which are listed in Table 1. 


Table I. RP-1 Reference Conditions 


mass flow (ibm/sec) 

64.29 

reference pressure (psi) 

37 

reference temperature (°R) 

530 

reference density (lmb/ft 3 ) 

50.3 

Reynolds number 

2.4 x 10 5 


Numerical Methodology 

The basic equations employed to describe the MPS 
feedline flow field are the multidimensional, general 
coordinate transport equations. A generalized form of 
these equations written in curvilinear coordinates is 
given by 

(l/J)(dpq/dt) = 3[-pU,q + |iG,j(3q/aj,)]/ dt, +(l/J)S q 

where J, Ui and Gy are written as: 

J = d(££j)/3(x,y) 

Ui = (Uj/J)(3£/dxj) 

Gy = (O^/Ox k ) (3£/dx k )/J, 

where q represents 1, u, v, w, h, k, and e, respectively. 
These are equations of continuity, x, y, and z 
momentum, enthalpy, turbulent kinetic energy, and 
turbulent kinetic energy dissipation rate. £ is the 
transformed curvilinear coordinate, p = (|i]+PtVcr H is the 
effective viscosity when the turbulent eddy viscosity 
concept is employed to model the turbulent flows, is the 
laminar viscosity; p 1 =pC M k 2 /e, is the turbulence eddy 
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viscosity and and cr H denote turbulence modeling 
constants. Turbulence modeling constants a 4 and source 
terms S H are given in Table 2. 

<J> is the energy dissipation function and P r represents the 
turbulent kinetic energy production term. C b C 2 , and C 3 
are model constants for the two-equation turbulence 
models. 


Table 2. and S 4 of the Transport Equations 


q 


S, 

l 

1.00 

0 

u 

1.00 

•p. + V[n( Uj)J -2/3(nVuj), 

v 

1.00 

-Py + V[n( Uj) y ] -2/3(|iVuj)j, 

w 

1.00 

-Pz + V[|4 a j) z ] -2/3(|iVu ,) z 

h 

0.95 

Dp/Dt + <)> 

k 

0.89 

P(P,-e) 

£ 

1.15 

p(e/k)(C,P r - C 2 e +CjP r 2 /e) 


To solve the system of nonlinear, coupled partial 
differential equations, finite difference approximations 
are used to establish a system of linearized algebraic 
equations. A relaxation solution procedure is then 
employed to couple the equations. An adaptive upwind 
scheme is employed to approximate the convective 
terms of the momentum, energy and continuity 
equations; the scheme is based on second- and fourth- 
order central differencing with artificial dissipation. 
First-order upwinding is used for species and turbulence 
equations since the parameters involved are positive 
quantities. Different eigenvalues are used for weighting 
the dissipation terms depending on the conserved 
quantity being evaluated in order to give correct 
diffusion fluxes near wall boundaries. Details of the 
present numerical methodology is given in [4, 5]. 

Discussion of Results 

The velocity profile just inside the pipe surface is shown 
in Figure 4. The flow separation can be seen only at the 
stagnation pocket to the right of the pressure- 
compensated bellow. Also, there is flow acceleration 
and deceleration through the elbows and neck-down 
regions. The close-up views of velocity vectors (Figure 
5) show the recirculation zones in the stagnation region, 
flow accelerations and decelerations over the bends, and 
uniform flow at the engine interface. A comparison of 
the velocity profiles along mutually perpendicular 
directions of the duct cross-section shows comparable 
velocity distortions at the engine interface for the new 
and the old designs (Figure 6). A more complete 
description of the flow field at the engine interface can 
be seen in the velocity contour plots for both designs 
(Figure 7). As seen in these plots, the new design 
indicates less flow variations at the engine interface. 


Figure 8 shows the static pressure on the pipe surface. 
These pressures are scaled to the reference pressure (37 
psi) at the engine interface. As expected, the bends 
generated pressure gradients normal to the bend radius, 

i.e. the pressure is low on the inner radius and higher on 
the outer radius of the elbow. The mass-averaged 
pressure distribution along the feedline model and its 
resulting pressure coefficients are shown in Figure 9 and 
Figure 10, respectively. Note that the elbows generated 
the secondary flows, which eventually contributed to the 
pressure losses, as seen in the pressure plots. The 
pressure drop between the inlet and the engine interface 
is about 4.6 psi for the new design, compared to 4. 1 psi 
in the previous design. However, this number may not 
be accurate since the model does not include the gimbal 
joints, fuel filters, and flow meters. These features are 
expected to have an impact on the pressure drop, but not 
much on the flow profiles. 

Conclusions 

In conclusion, the analyses indicate that the new duct 
design does not have a significant affect on the flow at 
the outlet of the feedline. It is also shown that the flow 
development is more favorable in the new design, 
although the pressure drop is slightly higher in this 
design. The introduction of the neck-down and the 
tighter bend did not grossly affect the flow and since the 
flow is fully developed at the engine interface, there is 
no need to use the flow straighteners at the engine inlet. 
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Figure 3: Old vs. New Design 
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Figure 4: Velocity Magnitude and Vectors for the RP-1 Fee dline 





X-34 RP-1 Feedline CFD Results 



relative distance across the pipe 
Figure 6: Velocity Profiles at the Engine Interface 


RP-1 Feedline CFD Results: Orbital’s Design RP-1 Feedline CFD Results (MSFCs Design) 

Non dimensional! zed Axial Velocity at the Engine Interface 

7 Nondimemionahzed Axial Velocity 



Figure 7: Engine Interface Velocity Contours 
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igure 8: 


RP 1 Feedline Mass Averaged Pressure 

Referenced to 37 psi. at the engine interface 


distance from the model inlet (inches) 

Figure 9: RP-1 Feedline Mass-Averaged Pressure Distribution 

q ret = 1.4915 psi 
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distance from the model inlet (inches) 

Figure 10: RP-1 Feedline Mass-Averaged Pressure Coefficient Distribution 








